Genetic introgression from commercial European pigs to the indigenous Chinese Lijiang breed and associated changes in phenotypes.

BACKGROUND: Gene flow is crucial for enhancing economic traits of livestock. In China, breeders have used hybridization strategies for decades to improve livestock performance. Here, we performed whole-genome sequencing of a native Chinese Lijiang pig (LJP) breed. By integrating previously published data, we explored the genetic structure and introgression of genetic components from commercial European pigs (EP) into the LJP, and examined the impact of this introgression on phenotypic traits. RESULTS: Our analysis revealed significant introgression of EP breeds into the LJP and other domestic pig breeds in China. Using a haplotype-based approach, we quantified introgression levels and compared EP to LJP and other Chinese domestic pigs. The results show that EP introgression is widely prevalent in Chinese domestic pigs, although there are significant differences between breeds. We propose that LJP could potentially act as a mediator for the transmission of EP haplotypes. We also examined the correlation between EP introgression and the number of thoracic vertebrae in LJP and identified VRTN and STUM as candidate genes for this trait. CONCLUSIONS: Our study provides evidence of introgressed European haplotypes in the LJP breed and describes the potential role of EP introgression on phenotypic changes of this indigenous breed.

AGIDB: a versatile database for genotype imputation and variant decoding across species.

The high cost of large-scale, high-coverage whole-genome sequencing has limited its application in genomics and genetics research. The common approach has been to impute whole-genome sequence variants obtained from a few individuals for a larger population of interest individually genotyped using SNP chip. An alternative involves low-coverage whole-genome sequencing (lcWGS) of all individuals in the larger population, followed by imputation to sequence resolution. To overcome limitations of processing lcWGS data and meeting specific genotype imputation requirements, we developed AGIDB (https://agidb.pro), a website comprising tools and database with an unprecedented sample size and comprehensive variant decoding for animals. AGIDB integrates whole-genome sequencing and chip data from 17 360 and 174 945 individuals, respectively, across 89 species to identify over one billion variants, totaling a massive 688.57 TB of processed data. AGIDB focuses on integrating multiple genotype imputation scenarios. It also provides user-friendly searching and data analysis modules that enable comprehensive annotation of genetic variants for specific populations. To meet a wide range of research requirements, AGIDB offers downloadable reference panels for each species in addition to its extensive dataset, variant decoding and utility tools. We hope that AGIDB will become a key foundational resource in genetics and breeding, providing robust support to researchers.

Genome-wide association study reveals genetic loci and candidate genes for meat quality traits in a four-way crossbred pig population.

Meat quality traits (MQTs) have gained more attention from breeders due to their increasing economic value in the commercial pig industry. In this genome-wide association study (GWAS), 223 four-way intercross pigs were genotyped using the specific-locus amplified fragment sequencing (SLAF-seq) and phenotyped for PH at 45 min post mortem (PH45), meat color score (MC), marbling score (MA), water loss rate (WL), drip loss (DL) in the longissimus muscle, and cooking loss (CL) in the psoas major muscle. A total of 227, 921 filtered single nucleotide polymorphisms (SNPs) evenly distributed across the entire genome were detected to perform GWAS. A total of 64 SNPs were identified for six meat quality traits using the mixed linear model (MLM), of which 24 SNPs were located in previously reported QTL regions. The phenotypic variation explained (PVE) by the significant SNPs was from 2.43% to 16.32%. The genomic heritability estimates based on SNP for six meat-quality traits were low to moderate (0.07-0.47) being the lowest for CL and the highest for DL. A total of 30 genes located within 10 kb upstream or downstream of these significant SNPs were found. Furthermore, several candidate genes for MQTs were detected, including pH45 (GRM8), MC (ANKRD6), MA (MACROD2 and ABCG1), WL (TMEM50A), CL (PIP4K2A) and DL (CDYL2, CHL1, ABCA4, ZAG and SLC1A2). This study provided substantial new evidence for several candidate genes to participate in different pork quality traits. The identification of these SNPs and candidate genes provided a basis for molecular marker-assisted breeding and improvement of pork quality traits.

Molecular genetic characterization and meat-use functional gene identification in Jianshui yellow-brown ducks through combined resequencing and transcriptome analysis.

The Jianshui yellow-brown duck is a unique country-specific waterfowl species in Yunnan Province, well known for its tender meat. However, there is a lack of comprehensive systematic research on the molecular genetic characteristics, especially germplasm resources and economic traits, of the Jianshui yellow-brown ducks. This study investigated the molecular genetic characteristics of Jianshui yellow-brown ducks, compared their selection signals with those of ancestral mallard and meat-type Pekin ducks, and identified genes specific to their meat-use performance. Furthermore, this study also evaluated the breeding potential for its meat performance. In this study, phylogenetic trees, PCA and Admixture analysis were used to investigate the population genetic structure among local duck breeds in China; population genetic differentiation index (Fst), nucleotide diversity and Tajima's D were used to detect selected loci and genes in the population of Jianshui yellow-brown ducks; and transcriptome technology was used to screen for differentially expressed genes in the liver, sebum and breast muscle tissues, and finally, the results of the genome selection signals and transcriptome data were integrated to excavate functional genes affecting the meat performance of the Jianshui yellow-brown ducks. The results of the genetic structure of the population showed that Jianshui yellow-brown ducks were clustered into a separate group. Selection signal analysis indicated significant selection pressure on certain genes related to meat characteristics (ELOVL2, ELOVL3, GDF10, VSTM2A, PHOSPHO1, and IGF2BP1) in both Jianshui yellow-brown ducks and mallards. Transcriptomic data analysis suggested that ELOVL3, PHOSPHO1, and GDF10 are vital candidate genes influencing meat production and quality in Jianshui yellow-brown ducks. A comparison of selection signals between Jianshui yellow-brown ducks and Pekin ducks revealed only 21 selected genes in the Jianshui yellow-brown duck population, and no significant genes were related to meat traits. Moreover, whole-genome resequencing data suggested that the Jianshui yellow-brown duck represents a unique category with distinct genetic mechanisms. Through selection signaling and transcriptomic approaches, we successfully screened and identified important candidate genes affecting meat traits in Jianshui yellow-brown ducks. Furthermore, the Jianshui yellow-brown duck has good potential for improved meat performance, highlighting the need for further improvement.

Genome-Wide Association Study of Growth Traits in a Four-Way Crossbred Pig Population.

Growth traits are crucial economic traits in the commercial pig industry and have a substantial impact on pig production. However, the genetic mechanism of growth traits is not very clear. In this study, we performed a genome-wide association study (GWAS) based on the specific-locus amplified fragment sequencing (SLAF-seq) to analyze ten growth traits on 223 four-way intercross pigs. A total of 227,921 highly consistent single nucleotide polymorphisms (SNPs) uniformly dispersed throughout the entire genome were used to conduct GWAS. A total of 53 SNPs were identified for ten growth traits using the mixed linear model (MLM), of which 18 SNPs were located in previously reported quantitative trait loci (QTL) regions. Two novel QTLs on SSC4 and SSC7 were related to average daily gain from 30 to 60 kg (ADG30-60) and body length (BL), respectively. Furthermore, 13 candidate genes (ATP5O, GHRHR, TRIM55, EIF2AK1, PLEKHA1, BRAP, COL11A2, HMGA1, NHLRC1, SGSM1, NFATC2, MAML1, and PSD3) were found to be associated with growth traits in pigs. The GWAS findings will enhance our comprehension of the genetic architecture of growth traits. We suggested that these detected SNPs and corresponding candidate genes might provide a biological foundation for improving the growth and production performance of pigs in swine breeding.

Genome-wide association study identifying genetic variants associated with carcass backfat thickness, lean percentage and fat percentage in a four-way crossbred pig population using SLAF-seq technology.

BACKGROUND: Carcass backfat thickness (BFT), carcass lean percentage (CLP) and carcass fat percentage (CFP) are important to the commercial pig industry. Nevertheless, the genetic architecture of BFT, CLP and CFP is still elusive. Here, we performed a genome-wide association study (GWAS) based on specific-locus amplified fragment sequencing (SLAF-seq) to analyze seven fatness-related traits, including five BFTs, CLP, and CFP on 223 four-way crossbred pigs. RESULTS: A total of 227, 921 highly consistent single nucleotide polymorphisms (SNPs) evenly distributed throughout the genome were used to perform GWAS. Using the mixed linear model (MLM), a total of 20 SNP loci significantly related to these traits were identified on ten Sus scrofa chromosomes (SSC), of which 10 SNPs were located in previously reported quantitative trait loci (QTL) regions. On SSC7, two SNPs (SSC7:29,503,670 and rs1112937671) for average backfat thickness (ABFT) exceeded 1% and 10% Bonferroni genome-wide significance levels, respectively. These two SNP loci were located within an intron region of the COL21A1 gene, which was a protein-coding gene that played an important role in the porcine backfat deposition by affecting extracellular matrix (ECM) remodeling. In addition, based on the other three significant SNPs on SSC7, five candidate genes, ZNF184, ZNF391, HMGA1, GRM4 and NUDT3 were proposed to influence BFT. On SSC9, two SNPs for backfat thickness at 6-7 ribs (67RBFT) and one SNP for CLP were in the same locus region (19 kb interval). These three SNPs were located in the PGM2L1 gene, which encoded a protein that played an indispensable role in glycogen metabolism, glycolysis and gluconeogenesis as a key enzyme. Finally, one significant SNP on SSC14 for CLP was located within the PLBD2 gene, which participated in the lipid catabolic process. CONCLUSIONS: A total of two regions on SSC7 and SSC9 and eight potential candidate genes were found for fatness-related traits in pigs. The results of this GWAS based on SLAF-seq will greatly advance our understanding of the genetic architecture of BFT, CLP, and CFP traits. These identified SNP loci and candidate genes might serve as a biological basis for improving the important fatness-related traits of pigs.

Genetic Diversity and Selection Signatures Within Diannan Small-Ear Pigs Revealed by Next-Generation Sequencing.

Genetic characterization of Chinese indigenous pig breeds is essential to promote scientific conservation and sustainable development of pigs. Here, we systematically surveyed the genomes of 75 unrelated Diannan small-ear (DSE) pigs from three diverse regions (Yingjiang County, Jinping County, and Sipsongpanna in Yunnan Province) to describe their population structures, genetic diversity, inbreeding coefficients, and selection signatures. First, these individuals were sequenced and genotyped using the genome reducing and sequencing (GGRS) protocol. A total of 438,038 autosomal single-nucleotide polymorphisms (SNPs) were obtained and used for subsequent statistical analysis. The results showed that these DSE pigs were clearly differentiated into three separate clades revealed by the population structure and principal component analysis, which is consistent with their geographical origins. Diannan small-ear pigs owned lower genetic diversity when compared with some other pig breeds, which demonstrated the need to strengthen the conservation strategies for DSE pigs. In addition, the inbreeding coefficients based on runs of homozygosity (ROH) length (F ROH) were calculated in each ROH length categories, respectively. And the results indicated that the ancient (up to 50 generations ago) inbreeding had greater impacts than recent (within the last five generations) inbreeding within DSE pigs. Some candidate selection signatures within the DSE pig population were detected through the ROH islands and integrated haplotype homozygosity score (iHS) methods. And genes associated with meat quality (COL15A1, RPL3L, and SLC9A3R2), body size (PALM2-AKAP2, NANS, TRAF7, and PACSIN1), adaptability (CLDN9 and E4F1), and appetite (GRM4) were identified. These findings can help to understand the genetic characteristics and provide insights into the molecular background of special phenotypes of DSE pigs to promote conservation and sustainability of the breed.

Bacterial magnetic particles-polyethylenimine vectors deliver target genes into multiple cell types with a high efficiency and low toxicity.

Bacterial magnetic particles (BMPs) are biosynthesized magnetic nano-scale materials with excellent dispersibility and biomembrane enclosure properties. In this study, we demonstrate that BMPs augment the ability of polyethylenimine (PEI) to deliver target DNA into difficult-to-transfect primary porcine liver cells, with transfection efficiency reaching over 30%. Compared with standard lipofection and polyfection, BMP-PEI gene vectors significantly enhanced the transfection efficiencies for the primary porcine liver cells and C2C12 mouse myoblast cell lines. To better understand the mechanism of magnetofection using BMP-PEI/DNA vectors, transmission electron microscopy (TEM) images of transfected Cos-7, HeLa, and HEP-G2 cells were observed. We found that the BMP-PEI/DNA complexes were trafficked into the cytoplasm and nucleus by way of vesicular transport and endocytosis. Our study builds support for the versatile BMP-PEI vector transfection system, which might be exploited to transfect a wide range of cell types or even to reach specific targets in the treatment of disease. KEY POINTS: • We constructed a BMP-PEI gene delivery vector by combining BMPs and PEI. • The vector significantly enhanced transfection efficiencies in eukaryotic cell lines. • The transfection mechanism of this vector was explained in our study.

Convergent genomic signatures of high-altitude adaptation among domestic mammals.

Abundant and diverse domestic mammals living on the Tibetan Plateau provide useful materials for investigating adaptive evolution and genetic convergence. Here, we used 327 genomes from horses, sheep, goats, cattle, pigs and dogs living at both high and low altitudes, including 73 genomes generated for this study, to disentangle the genetic mechanisms underlying local adaptation of domestic mammals. Although molecular convergence is comparatively rare at the DNA sequence level, we found convergent signature of positive selection at the gene level, particularly the EPAS1 gene in these Tibetan domestic mammals. We also reported a potential function in response to hypoxia for the gene C10orf67, which underwent positive selection in three of the domestic mammals. Our data provide an insight into adaptive evolution of high-altitude domestic mammals, and should facilitate the search for additional novel genes involved in the hypoxia response pathway.

Natural selection on TMPRSS6 associated with the blunted erythropoiesis and improved blood viscosity in Tibetan pigs.

Tibetan pigs, indigenous to Tibetan plateau, are well adapted to hypoxia. So far, there have been not any definitively described genes and functional sites responsible for hypoxia adaptation for the Tibetan pig. The whole genome-wide association studies in human suggested that genetic variations in TMPRSS6 was associated with hemoglobin concentration (HGB) and red cell counts (RBC). Here we conducted resequencing of the nearly entire genomic region (40.1 kb) of the candidate gene TMPRSS6 in 40 domestic pigs and 40 wild boars along continuous altitudes and identified 708 SNPs, in addition to an indel (CGTG/----) in the intron 10. We conduct the CGTG indel in 838 domestic pigs, both the CGTG deletion frequency and the pairwise r2 linkage disequilibrium showed an increase with elevated altitudes, suggesting that TMPRSS6 has been under Darwinian positive selection. As the conserved core sequence of hypoxia-response elements (HREs), the deletion of CGTG in Tibetan pigs decreased the expression levels of TMPRSS6 mRNA and protein in the liver revealed by real-time quantitative PCR and western blot, respectively. We compared domestic pigs and Tibetan pigs living continuous altitudes, found that the blood-related traits with the increase of altitude, however, the HGB did not increase with the elevation in Tibetan pigs. Genotype association analysis results dissected a genetic effect on reducing HGB by 13.25 g/L in Gongbo'gyamda Tibetan pigs, decreasing mean corpuscular volume (MCV) by 4.79 fl in Diqing Tibetan pigs. In conclusion, the CGTG deletion of TMPRSS6 resulted in lower HGB and smaller MCV, which could reflect a blunting erythropoiesis and improving blood viscosity as well as erythrocyte deformability. It remains to be determined whether a blunting of erythropoiesis for TMPRSS6 or others genetic effects are the physiological adaptations among Tibetan pigs.

Haplotype diversity in mitochondrial DNA reveals the multiple origins of Tibetan horse.

The Tibetan horse is a species endemic to the Tibetan plateau, with considerable economic value in the region. However, we currently have little genetic evidence to verify whether the breed originated in Tibet or if it entered the area via an ancient migratory route. In the present study, we analyzed the hypervariable segment I sequences of mitochondrial DNA (mtDNA) in 2,050 horses, including 290 individuals from five Tibetan populations and 1,760 from other areas across Asia. Network analysis revealed multiple maternal lineages in the Tibetan horse. Component analysis of sub-lineage F3 indicated that it decreased in frequency from east to west, a trend reflected both southward and northward from Inner Mongolia. Analysis of population genetics showed that the Deqen horse of eastern Tibet was more closely related to the Ningqiang horse of northern China than to other Tibetan horses or the Yunnan horse. These results indicated that the Tibetan horse migrated first from Central Asia to Mongolia, moved south to eastern Tibet (near Deqen), then finally westward to other regions of Tibet. We also identified a novel lineage K that mainly comprises Tibetan and Yunnan horses, suggesting autochthonous domesticated origin for some Tibetan horse breeds from local wild horses. In conclusion, our study demonstrated that modern Tibetan horse breeds originated from the introgression of local wild horses with exotic domesticated populations outside China.

Attaching Biosynthesized Bacterial Magnetic Particles to Polyethylenimine Enhances Gene Delivery Into Mammalian Cells.

Gene transfection using bacterial magnetic particles (BMPs)-polyethylenimine (PEI) has become increasingly prevalent; however, relatively little effort has been made to optimize the protocol for preparing these complexes with the aim of improving their transfection efficiency. Here, we report a procedure for constructing BMPs-PEI/DNA complexes that results in improved transfection efficiency, reduced cytotoxicity and shorter procedure times for both complex formation and transfection over current methods. BMPs-PEI/DNA complexes mixed using ultrasonication yielded beads that were 10.2% more efficient at transfecting HeLa cells than complexes made by mechanical vortexing. Phosphate-buffered saline (PBS) proved to be a superior solvent for BMPs-PEI/DNA and PEI/DNA complexes, and the transfection efficiencies in HeLa cells were 54.75% and 46.01%, respectively. Comparable levels of transfection were achieved after 10 min of incubation with low-dose BMPs-PEI/DNA complexes versus 4 h with standard PEI/DNA complexes. BMPs stored in PBS have an average transfection efficiency that is 5% greater than those stored in physiological salt solutions. Cell morphology and cytotoxicity analyses demonstrated that the biosynthesized BMPs lessened the cytotoxicity of PEI to cells. Our results provide an optimized protocol for BMPs-PEI/DNA complex construction and gene transfer in vitro.

Whole-genome sequencing of six dog breeds from continuous altitudes reveals adaptation to high-altitude hypoxia.

The hypoxic environment imposes severe selective pressure on species living at high altitude. To understand the genetic bases of adaptation to high altitude in dogs, we performed whole-genome sequencing of 60 dogs including five breeds living at continuous altitudes along the Tibetan Plateau from 800 to 5100 m as well as one European breed. More than 150× sequencing coverage for each breed provides us with a comprehensive assessment of the genetic polymorphisms of the dogs, including Tibetan Mastiffs. Comparison of the breeds from different altitudes reveals strong signals of population differentiation at the locus of hypoxia-related genes including endothelial Per-Arnt-Sim (PAS) domain protein 1 (EPAS1) and beta hemoglobin cluster. Notably, four novel nonsynonymous mutations specific to high-altitude dogs are identified at EPAS1, one of which occurred at a quite conserved site in the PAS domain. The association testing between EPAS1 genotypes and blood-related phenotypes on additional high-altitude dogs reveals that the homozygous mutation is associated with decreased blood flow resistance, which may help to improve hemorheologic fitness. Interestingly, EPAS1 was also identified as a selective target in Tibetan highlanders, though no amino acid changes were found. Thus, our results not only indicate parallel evolution of humans and dogs in adaptation to high-altitude hypoxia, but also provide a new opportunity to study the role of EPAS1 in the adaptive processes.

Investigation on the transcription factors of porcine 3ß-hydroxysteroid dehydrogenase and 17ß-hydroxysteroid dehydrogenase genes.

The enzymes 3ß-hydroxysteroid dehydrogenase (3ßHSD) and 17ß-hydroxysteroid dehydrogenase (17ßHSD) regulate the steroid metabolism in mammals. In this study, we aimed to characterize the steroid related transcription factors at the 5' flanking region of these two genes. A series of 5' deletions of approximately 1 kb of 5'-flanking region on both genes were fused to a pGL3 basic vector containing firefly luciferase cDNA, and then transfected to human hepatocellular liver carcinoma cell line (HepG2). Luciferase activity assay indicated the region from -574 to -617 bp of the 3ßHSD1 promoter, and from -850 to -868 bp of 17ßHSD7 promoter induced the highest luciferase activity. A putative transcription factor, i.e. the proline and acidic amino acid-rich basic leucine zipper (PAR/bZIP) family of 3ßHSD1 gene, and three-amino acid loop extension (TALE) homeodomain class of 17ßHSD7 were identified respectively by sequence homology. Gel shift assay further confirmed the binding capacity of the putative elements to nuclear extract. Our study gives new insights to the transcriptional regulation of 3ßHSD1 and 17ßHSD7 and further hints to their involvement in steroid metabolism.

Investigation of testosterone, androstenone, and estradiol metabolism in HepG2 cells and primary culture pig hepatocytes and their effects on 17ßHSD7 gene expression.

Steroid metabolism is important in various species. The accumulation of androgen metabolite, androstenone, in pig adipose tissue is negatively associated with pork flavor, odour and makes the meat unfit for human consumption. The 17ß-hydroxysteroid dehydrogenase type 7 (17ßHSD7) expressed abundantly in porcine liver, and it was previously suggested to be associated with androstenone levels. Understanding the enzymes and metabolic pathways responsible for androstenone as well as other steroids metabolism is important for improving the meat quality. At the same time, metabolism of steroids is known to be species- and tissue-specific. Therefore it is important to investigate between-species variations in the hepatic steroid metabolism and to elucidate the role of 17ßHSD7 in this process. Here we used an effective methodological approach, liquid chromatography coupled with mass spectrometry, to investigate species-specific metabolism of androstenone, testosterone and beta-estradiol in HepG2 cell line, and pig cultured hepatocytes. Species- and concentration-depended effect of steroids on 17ßHSD7 gene expression was also investigated. It was demonstrated that the investigated steroids can regulate the 17ßHSD7 gene expression in HepG2 and primary cultured porcine hepatocytes in a concentration-dependent and species-dependent pattern. Investigation of steroid metabolites demonstrated that androstenone formed a 3'-hydroxy compound 3ß-hydroxy-5α-androst-16-ene. Testosterone was metabolized to 4-androstene-3,17-dione. Estrone was found as the metabolite for ß-estradiol. Inhibition study with 17ßHSD inhibitor apigenin showed that apigenin didn't affect androstenone metabolism. Apigenin at high concentration (50 µM) tends to inhibit testosterone metabolism but this inhibition effect was negligible. Beta-estradiol metabolism was notably inhibited with apigenin at high concentration. The study also established that the level of testosterone and ß-estradiol metabolites was markedly increased after co-incubation with high concentration of apigenin. This study established that 17ßHSD7 is not the key enzyme responsible for androstenone and testosterone metabolism in porcine liver cells.